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Abstract 

In this paper, we revisit the claim that the Eulerian and quasi-Lagrangian same time 
correlation tensors are equal. This statement allows us to transform the results of 
an MSR quasi-Lagrangian statistical theory of hydrodynamic turbulence back to 
the Eulerian representation. We define a hierarchy of homogeneity symmetries be- 
tween incremental homogeneity and global homogeneity. It is shown that both the 
elimination of the sweeping interactions and the derivation of the 4/5-law require a 
homogeneity assumption stronger than incremental homogeneity but weaker than 
global homogeneity. The quasi-Lagrangian transformation, on the other hand, re- 
quires an even stronger homogeneity assumption which is many-time rather than 
one-time but still weaker than many-time global homogeneity. We argue that it is 
possible to relax this stronger assumption and still preserve the conclusions derived 
from theoretical work based on the quasi-Lagrangian transformation. 

Key words: Turbulence, local homogeneity, quasi-Lagrangian, sweeping 
interactions 

PACS: 47.27.Ak, 47.27.Gs, 47.27.Jv 



1 Introduction 



A remarkable feature of hydrodynamic turbulence in three dimensions is that it exhibits uni- 
versal self-similarity properties at small length scales independently of the forcing mechanism that 
operates at larger length scales. The self-similar nature of turbulence was noticed by Richardson 
[1] who suggested that large vortices will generate increasingly smaller vortices until they become 
hydrodynamically stable and then get dissipated by viscosity. Kolmogorov [2, 3] conjectured that 
for length scales r between the forcing scale Iq and the dissipation scale ry, the structure functions 
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Sni^,re) will be independent of Iq and rj, and, as was pointed out by Batchelor [4], this conjecture 
implies that Sn{x,re) satisfy the following power laws: 



5„(x, re) = ({[u(x + re, t) - u(x, t)] ■ e}") = C„(£r)"/^ 



(1) 



Here e is a unit vector, and e equals the rate of energy injection into the fluid, the energy flux in 
the cascade of energy from large scales to small scales, and the rate of energy dissipation at small 
scales. The constant C„ was believed to be universal, but in fact it is not (except for n = 3) and it is 
dependent on the forcing spectrum. Prom the above, the energy spectrum E(k) for i^^ ^ A; <C ri~^ 
can be shown to satisfy 



This prediction was confirmed for the first time in 1962 [5,6], and today, with modern computers, 
it is routinely reproduced in numerical simulations. It has since come to light [7, 8] that there 
exist departures from Kolmogorov scaling laws for the higher order structure functions (known as 
intermittency corrections), and Kolmogorov (with Oboukhov) [9, 10] was in fact the first to propose 
revisions of his original theory. The correct expression for Sn{r) has the form 



where d are scaling exponents to be determined. The challenge here has been to develop theoretical 
understanding that can account for this energy cascade with a logical argument that begins from 
the underlying governing Navier-Stokes equations. It is not only a matter of calculating the scaling 
exponents Cn- The robustness of the scaling of the energy spectrum needs to be explained, and the 
universality of the scaling exponents themselves is in fact still an open question. 

The energy cascade from large scales to small scales is driven by the nonlinear term of the 
Navier-Stokes equations, and it is often explained as an effect of the vortex stretching and tilting 
caused by that term. However, the same term is also responsible for a sweeping interaction whereby 
a vortex is swept altogether from one location to another with minimal distortion. Implicit in the 
idea of an energy cascade is the assumption that these sweeping interactions have a negligible effect 
on the structure functions in the inertial range. It has therefore been necessary to use theoretical 
schemes that "eliminate" sweeping [11]. The goal of this paper is to call attention to the fact 
that these schemes do not prove that sweeping is negligible; they only introduce the assumption 
that it is so. Recent doubts [12] concerning the consistency of the local homogeneity framework 
are directly linked with this problem of rigorously eliminating the sweeping interactions, and also 
with the problem of formulating a reasonable definition of local homogeneity. We will also make a 
conjecture, and explore its plausibility, which, if shown to be true, would establish the assumption 
that the sweeping interactions are negligible in the inertial range on a firmer ground. 

It should be noted that a strictly rigorous mathematical theory based exclusively on the Navier- 
Stokes equations is a very difficult task. For this reason, it is necessary to tolerate unproven as- 
sumptions as hypotheses, as long as such assumptions can be reasonably supported by physical 
arguments, or by experiment. It is within a specific framework of reasonable assumptions, which 
will be defined in a moment, that we claim that sweeping elimination procedures still do not prove 
that sweeping interactions are negligible. 

The argument of this paper, summarily, is the following. First, we show that the elimination of 
the sweeping interactions as well as the derivation of the 4/5-law requires a homogeneity assumption 
stronger than the assumption of incremental homogeneity, as envisioned by Frisch [7]. Second, 
we show that using the quasi-Lagrangian formulation of Belinicher and L'vov to eliminate the 



E{k) = C£2/3A;-V3. 



(2) 



Sn{r) = Cn{erT'\r/iof--^'\ 



(3) 
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sweeping interactions requires an even stronger homogeneity assumption which involves many-time 
correlations instead of one-time correlations. We conclude with a discussion of the implications of 
this argument on the utility of the quasi-Lagrangian formulation. Specifically, we will show that 
despite this apparent shortcoming of the quasi-Lagrangian formulation, the theoretical work based 
on it can still be used as a foundation for a physically useful theory, along the lines of the Prisch 
framework, provided that certain considerations are taken into account. Furthermore, incremental 
homogeneity is in fact a consistent framework, provided that the sweeping interactions can be 
eliminated in a more rigorous manner. 

The paper is organized as follows. In section 2, we review the theoretical developments that gave 
rise to the issue of the sweeping interactions, and discuss the assumptions underlying most efforts to 
understand the energy cascade from a theoretical point of view. The 4/5- law is discussed in section 3 
and the quasi-Lagrangian formulation in section 4. The implications of our argument for the theories 
that use the quasi-Lagrangian formulation as a foundation are discussed in section 5, and the paper 
is concluded in section 6. Appendix A reviews how the quasi-Lagrangian formulation eliminates 
the sweeping interactions. Appendix B presents a more complete account of the calculation of a 
functional determinant originally given by L'vov and Procaccia [13]. In appendix C we evaluate the 
contribution of the sweeping interactions in closed form for the case of passive random gaussian 
delta-time-decorrelated sweeping. 

2 Theoretical background 

We begin with reviewing the development of the ideas that form the theoretical foundation of 
certain recent attempts to understand the universal behavior of turbulence. The problem of the 
sweeping interactions and its resolution is an essential part of this theoretical foundation. Then we 
discuss the set of assumptions that are widely accepted on physical grounds. Because the argument 
of this paper requires simultaneous consideration of a wide range of interdependent topical areas, 
this overview will help by providing the reader the broader context against which the argument and 
its implications on the theoretical foundations of turbulence will be discussed later in our paper. 
This overview represents strictly my personal philosophical point of view. A more comprehensive 
and unbiased review of theoretical three-dimensional turbulence is already available in the literature 
[7,11,14, 15]. 

2.1 Theoretical approaches to turbulence 

The foundation on which recent successful theoretical work was accomplished on the prob- 
lem of the direct energy cascade rests on the following essential ideas: The first critical idea is 
the framework of globally homogeneous and isotropic turbulence introduced by Taylor [16-19] and 
popularized by Batchelor [20] . Within that framework there have been numerous attempts to model 
turbulence using closure models [21]. The second critical idea, due to Kraichnan, is his discovery 
that such models are not realizable because they predict negative values for the energy spectrum 
[22]. Kraichnan counterproposed a different closure model [23,24], the direct interaction approxi- 
mation (DIA), with the unique feature that it makes use of response functions. Disagreement with 
experimental predictions prompted Kraichnan to call attention to the problem of sweeping inter- 
actions [25], and to revise his earlier model. The new model [26], the Lagrangian history direct 
interaction approximation (LHDIA), was one of the first models to make predictions in agreement 
with experiment [27]. A review of Kraichnan's work was given by Leslie [21]. It is fair to say that 
LHDIA was the first successful theory of three-dimensional turbulence. Unfortunately, it was not 
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clear how to generalize LHDIA, which was a first order approximation, to higher orders, and as a 
result, further development of this theoretical program was not possible for many years. 

Parallel to these efforts, there have also been attempts to construct exact mathematical theories 
of turbulence based on functional calculus. The first such formulation was given by Hopf [28], and 
an equivalent reformulation in terms of path integrals by Rosen [29, 30]. Novikov [31] modified the 
Hopf formalism to include a gaussian delta correlated stochastic forcing, intended to model the 
hydrodynamic instability responsible for turbulence. An interesting application of this formalism 
is the more rigorous and powerful reformulations of the original dimensional analysis arguments 
used by Kolmogorov [32,33]. Its main disadvantage is that it restricts the statistical description 
to one-time velocity correlations. A generalization to include many-time velocity correlations was 
given by Lewis and Kraichnan [34]; however even that is inadequate because it does not include 
response functions. 

The essential idea of the definitive approach was introduced by Wyld [35]. The main result 
is that Feynman diagrams can be used to generalize DIA to higher orders, and that DIA itself 
is essentially a one-loop line-renormalized diagrammatic theory. A generalization of this scheme 
to a wider range of dynamical systems was given by Martin, Siggia, and Rose [36], although, as 
they themselves explained, without a sufficiently rigorous justification. Phythian [37] used Feynman 
path integrals to reformulate the MSR theory, and showed that it can be justified for dynamical 
systems that are local in time and first-order in time. An assumption implicit in this argument is 
that the dynamical system has a unique solution for all time. This claim has not been proven for 
the Navier-Stokes equations in three dimensions, however global regularity, as a matter of fact, can 
be proved rigorously [38, 39] if the diffusion term in the Navier-Stokes equations is replaced with a 
hyperdiffusion term like uV^Ua- A pedagogical introduction to MSR theory was given recently by 
L'vov and Procaccia [40] and Eyink [41], and a careful review of the mathematical foundations of 
the theory itself is given in the paper by Andersen [42] (also see references therein). 

Unfortunately, the MSR formalism could not be applied to generalize Kraichnan's more suc- 
cessful LHDIA theory because the Navier-Stokes cr^uations in the Lagrangian representation are not 
local in time. Eventually, a way was discovered around this difficulty, thus breaking the deadlock 
that has been plagueing theory for decades. It involves combining the MSR formalism with renor- 
malization schemes that eliminate the sweeping interactions. The first such scheme was introduced 
by Yakhot [43], and another by Belinicher and LVov [44,45]. Combined with the MSR formalism, 
one has a rather solid foundation for further theoretical work. It is these schemes, and the nature 
of the assumptions that they implicitly introduce, that will concern us in this paper. 

Since 1995, there have been some very remarkable developments in this direction: L'vov and 
Procaccia have used the quasi-lagrangian renormalization scheme [44, 45] to formulate a diagram- 
matic theory [13, 46, 47] that generalized Kraichnan's DIA to all orders. It was shown that as long 
as the theory is truncated to finite order, it predicts agreement with Kolmogorov's theory and the 
absence of intermittency corrections [13]. It was also shown that if the theory is not truncated, there 
is a critical divergence that does lead to intermittency corrections [46]. L'vov and Procaccia et al 
also formulated a nonperturbative theory [48-51] based on the fusion rules which are predicted 
by the underlying diagrammatic theory. This theory has been used to derive a nonperturbative 
method [52-54] and a perturbative method [55] for calculating the scaling exponents (^n- The per- 
turbative method [55] has been used successfully to calculate Cn for all n accessible to experimental 
measurement, but it requires that the deviation of C2 from the Kolmogorov prediction 2/3, which 
is the small parameter, be already known. This ability of the L'vov-Procaccia theory to predict the 
existence of intermittency is a significant accomplishment. A partial review of these developments 
was given in [11]. The non-perturbative theory has also led to a clearer understanding of local 
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isotropy [56, 57]. 

It is worth mentioning that there exists an entirely different theoretical approach to the problem 
based on renormalization group methods. A detailed review is given in [7, 58, 59] and some relevant 
criticism in Refs. [60, 61). There are two interesting points of convergence between renormalization 
group methods and the theories reviewed previously. First, Eyink [62, 63] employed the renormal- 
ization group method to derive the fusion rules under certain assumptions both for shell models of 
turbulence and for hydrodynamic turbulence itself. The fusion rules are a crucial element in both 
the perturbative and the non-perturbative theories of L'vov and Procaccia. Second, Giles [64] used 
the renormalization group method to calculate the scaling exponents Cn, without relying on any 
experimental input, contrary to the paper [55]. In this calculation, the sweeping interactions were 
eliminated using the scheme by Yakhot [43]. A comparative study of the two approaches would 
help further progress. 

2.2 The hypotheses that underlie MSR theory 

In all the theoretical work that has been reviewed above, it is assumed that the Navier-Stokes 
equations have a unique solution, that there exists hydrodynamic instability leading to turbulence, 
and that this instability can be modeled with stochastic forcing acting at large scales. These assump- 
tions are introduced implicitly simply by employing the MSR formalism. Although they are widely 
accepted on physical grounds, there has also been substantial effort to deal with them rigorously. 

An overview of the mathematical results on the existence and uniqueness of solutions to the 
Navier-Stokes equations is given in ref. [65, 66] and references therein. Briefly, in two dimensions 
the existence and uniqueness of strong solutions has been shown rigorously. In three dimensions it 
has been shown that weak solutions exist, but not that they are unique. It has also been shown that 
if strong solutions exist, they will have to be unique, but it has not been shown that such strong 
solutions do in fact exist. The underlying physical issue is whether the velocity field will develop 
singularities by vortex stretching as it is evolved by the Navier-Stokes equations. 

It is fortunate that this issue does not arise in numerical simulations because the finiteness of 
the resolution prevents singularities from developing. As long as the smallest resolved length scale is 
smaller by order of magnitudes than the Kolmogorov microscale, the finite resolution approximation 
of the Navier-Stokes equations models hydrodynamic turbulence quite adequately. Furthermore, the 
energy cascade, which is very robust, will not allow any of the Fourier modes to blow out, since 
all the incoming energy will be transferred to the dissipation range, where it will be disposed of 
efficiently, given adequate numerical resolution. Another benefit of the finite resolution model is 
that the path integrals of the corresponding MSR theory are mathematically rigorous. 

It should be noted that the Navier-Stokes equations themselves are not obviously more realistic 
than the finite resolution model because a "finite resolution" is imposed on fluid dynamics by Nature 
at the point where the existence of discrete molecules is important. Thus, if one introduces the 
assumption that the finite resolution approximation of the Navier-Stokes equations is a satisfactory 
physical model all by itself, then one may disregard the mathematical issues associated with the 
existence, uniqueness, and regularity of the solutions of the Navier-Stokes equations. This is not 
an unreasonable assumption in the inertial range of three-dimensional turbulence. We are on less 
solid ground with respect to the robustness of the cascades of two-dimensional turbulence, but the 
underlying mathematical issues do not arise in two dimensions. We do not wish to underestimate 
the importance of the mathematical issues of existence and uniqueness that remain open for 3D 
Navier-Stokes; we merely want to highlight the implicit assumption that one makes when one 
sidesteps these issues, as is done by every theory published to date. 
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Another very important issue which is "hidden under the rug" is proving tlie existence of 
turbulence itself as a consequence of the Navier-Stokcs equations. Unfortunately, the theoretical 
framework prescribed by the MSR theory cannot account, even in principle, for the existence 
of the hydrodynamic instability that causes turbulence. In the MSR framework, it is implicitly 
assumed that the effect of hydrodynamic instability can be modeled by a stochastic forcing term. 
The assumption can be justified if one demonstrates that the resulting stochastic behavior of the 
velocity field in the inertial range is invariant with respect to large-scale perturbations to the 
statistics of the forcing term. 

There is in fact an extension of MSR theory in terms of a supersymmetric path integral that 
includes two additional fermionic ghost fields [67-69]. The surprising result is that correlations 
involving these additional fields are related to the Lyapunov exponents [70] that quantify hydro- 
dynamic instability. It is therefore possible, in principle, to obtain statistical predictions from this 
framework with a deterministic forcing as input [71]. Whether this is in fact a practical approach 
remains to be seen. 

The assumptions described so far are needed to bring in the machinery of the MSR formalism. 
In order to employ the formalism to explain the universality of the direct energy cascade and 
calculate the intermittency corrections, it is necessary to hypothesize a mathematical description 
of the energy cascade and use that to narrow down the specific solution which is self-consistent. 
Prisch [7, 72] proposed a set of hypothesis consisting of assumptions of statistical symmetry (such as 
homogeneity, isotropy, self-similarity) and the additional assumption of anomalous dissipation, as an 
appropriate refinement of Kolmogorov's theory. The nature of the theoretical argument is to show 
that there is only a unique solution that can be admitted that satisfies the hypothesized statistical 
symmetries. A critical review of the assumed statistical symmetries, and local vs. incremental 
homogeneity in particular, is part of what concerns us in this paper. 

To summarize, we accept the following assumptions on physical or experimental grounds: first, 
there exists a unique solution to the Navier-Stokes equations that develops hydrodynamic instability 
for large Reynolds numbers; second, in the limit of fully developed turbulence, incremental homo- 
geneity and incremental isotropy (as defined by Prisch [7,72], and see section 3.2) are reinstated 
statistically, even if only asymptotically, for the velocity field; third, we accept the hypothesis that 
there exists an anomalous energy sink at small scales. These assumptions are a reasonable starting 
point for analytical theories of turbulence in three dimensions. 

3 Homogeneity and sweeping interactions 

The background on homogeneity is as follows: Taylor, Batchelor, Kraichnan, and others, have 
been willing to tolerate the assumption that turbulence is globally homogeneous and isotropic. 
However, it was suggested by Kolmogorov himself [2] that a far more realistic approach is to assume 
local homogeneity and local isotropy. Both frameworks have been reviewed by Monin and Yaglom 
[73]. Kolmogorov also emphasized the importance of studying stationary turbulence, corresponding 
to the forced-dissipative case, instead of the free decaying case. 

In recent work, Frisch [7, 72] proposed that Kolmogorov's second paper [3] leads to a refor- 
mulation of his theory along three assumptions: first, the assumption of local homogeneity and 
local isotropy (defined differently than by Kolmogorov, see section 3.2); second, an assumption of 
self-similarity; third, the assumption of an anomalous energy sink. Using the first and third assump- 
tion, according to Frisch, one derives the 4/5 law from which we obtain C3 = 1. From the second 
assumption we have (^n = nh. Combined, we obtain the prediction C,n = n/^. The assumption 
of self-similarity, used by Prisch, axiomatically excludes intermittency corrections to the scaling 
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exponents Consequently, the theoretical efforts to calculate the scaling exponents from "first 
principles" essentially aim to weaken this assumption while tolerating the other two assumptions. 

Some faith in the assumption of an anomalous energy sink, in particular, is based on recent 
evidence from numerical simulations [74] and theoretical evidence from the fusion rules [48, 50]. The 
assumption of local isotropy can be understood from the principle of linear superposition of the 
isotropic and anisotropic sectors of the symmetry group SO (3) [56,57]. Finally, the assumption of 
self-similarity can be understood via Z{h) covariance of the statistical theory [52-54]. This leaves 
then the assumption of local homogeneity. 



3.1 Hierarchical definitions of homogeneity 

Let Ua{^,t) be the Eulerian velocity field, and introduce the Eulerian velocity differences Wa- 
Wa{^,x' ,t) = Ua{x,t) - Ua (x' , t) . (4) 

The Eulerian generalized structure function is defined as the ensemble average of the product of 
such velocity differences 



F„"i"--'^({x,x'}„,t) 



Y[ W'afc(xfe,x'fe,t) 



.k=l 



(5) 



where {x, x'}„ is shorthand for a list of n position vectors. 

Originally, Frisch [7, 72] wrote his definitions of local homogeneity, local isotropy, and local 
stationarity using an "equivalence in law" relation. It should be noted that one should distinguish 
between many-time equivalence, that extends to many-time correlations, and one-time equivalence 
that applies only to one-time correlations. The clearest way to bring out this distinction is by 
defining the equivalence relation in terms of characteristic functionals defined as 



(6) 
(7) 



■^w'^'IP'*] = ^exp (^i J dxj dx' u;«(x,x',t)pa(x,x')) 
■^w''''*[p] = (|exp (^i J dxj dx' J i(;a(x,x',t)pa(x,x',i)) 



The structure functions can be evaluated from the characteristic functional by variational differen- 
tiation and setting p = 0. For example. 



i^r"'"""({x,x'}„,t) = 



S 



n 

TT-- 



[P,t] 



p=0 



(8) 



The difference between Z^'^ [p] and Z^'^ '* [p] , is that Z^'^ [p] contains information only about one- 
time correlations, whereas Z^'^ '*[p] contains information about many-time correlations as well. This 
is exploited to distinguish between many-time equivalence and one-time eqiiivalence. 
Definition 1 Consider two stochastic fields Ua(x,x',i) and it;a(x, x', i). The "equivalence in law" 
relations are defined as 



Va(x,x',t) w;a(x,x',t) <^ Z^''''[p,t\ = Z^''''[p,t\ Vp analytic 

p] Vp analytic. 



Va (X, x' , t) ~ '* Wa (x, x' , t) 



(9) 
(10) 
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Here, ^~ 



x,x ,t 



represents one-time equivalence, and '~' represents many-time equivalence. Thus, we 
can distinguish between one-time global homogeneity u G and many-time global homogeneity 
u G 



u G ?{ 
u G ?{* 



Ua{yi,t)^Ua{x + y,t), Vy GM'^ 
ti„(x,t)^*Ua(x + y,t), VyGM''. 



(11) 
(12) 



A detailed review of previous definitions of local homogeneity has been given by Hill [75]. To 
discuss local homogeneity more carefully, we introduce the following definitions: 

Definition 2 The velocity field u, as a stochastic field, is a member of the homogeneity class 
!Km{A) where yi C M*^ a region in W^, if and only if the ensemble average defined as 



.1=1 



.k=l 



(13) 



is invariant with respect to a space shift of its arguments xi,yk,y'k f^T all n > in the domain A, 



I.e. 



+ IZ(^A.yfc + ^A.y'fc) ^rn,n = 0, Vx/, y^, y'^ G A 



(14) 



fc=l 



Definition 3 The velocity field u is a member of the homogeneity class 'K^{A) where A C.W^ a 
region in W^, if and only if the ensemble average defined as 



p* = 

m.n — 



JJua(x;,tO 



.1=1 



Ylwp^{yk,y'k,t) 



.k=l 



(15) 



is invariant with respect to a space shift of its arguments xi,yk,y'k foi^ all n > in the domain A, 



I.e. 



dau^i + J^(%,yfc + 5/3fe,y' J Cn = 0, Vx^, y^, y';^ G A 



(16) 



k=l 



We also write "Km = 'Km{^'^) and "K*^ = "K^^iW^). The distinction between "KmiA) and "K^^iA) is 
that the former requires translational invariance on the one-time correlation tensor Fm,n, whereas 
the latter requires translational invariance on the many-time correlation tensor F^^, both over the 
domain A. 

We also define the following transfinite homogeneity classes: 



^U^) = n ^k{^) and %*JA) = fl :KI{A). 
fc€N fceN 



(17) 



In these homogeneity classes the ensemble average of any product of velocities multiplied with 
any product of velocity differences will be invariant under spatial shifting. Note that even this 
homogeneity class is weaker than global homogeneity. We will also distinguish between one-time 
global homogeneity u G J{ and many-time global homogeneity u G which are defined as 



u G 5{ <;=^ 'Uq(x, t) ~ Ua{x + y, t), Vy G 
ueO<* Ua{y:, t) ~ Uaipc + y, t), Vy G W^. 



(18) 
(19) 
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Remark 1 An immediate consequence of these definitions is that the homogeneity classes are hi- 
erarchically ordered, according to the following relations 



jf* c :k*ja) c :ki{a), yk g n, 

^a{A) C :Kb{A) A :KI{A) C :KI{A), Va, 6 G N : a > 6, 
^*a{A) C J{„(A), Va G N. 



(20) 

(21) 
(22) 
(23) 



3.2 Remarks on Kolmogorov's and Frisch's definition of local homogeneity 

The term "local homogeneity" is usually identified with the definition that was given by Kol- 
mogorov [2]. However, in his reformulation of the Kolmogorov 1941 theory, Frisch [7, 72] identified 
local homogeneity, local isotropy, and local stationarity with incremental homogeneity, incremental 
isotropy, and incremental stationarity. The definitions that he gave read: 

Locally stationary: Wai'x, x', i) Wa(x, x ,t + At) , VAt G M. 
Locally homogeneous: Waix, x', t) w-'q(x + y, x' + y, t) , Vy G M.'^. 

Locally isotropic: w^a(x, x', t) Wa(xo + ^(x — xq), xq + A(x' — xq), t) , \fA G SO{d). 

Using our notation, the condition of incremental homogeneity can be written as u G ^o{A). It 
should be stressed that Frisch postulated that these symmetries are valid asymptotically for space 
shifts and time shifts up to a relevant order of magnitude and proposed them as reasonable hy- 
potheses to be used as the basis for a modern reformulation of Kolmogorov's 1941 theory [7, 72]. 

To motivate his hypotheses, Frisch argues that homogeneity, isotropy, and time invariance are 
satisfied by the Navier-Stokes equations and they are violated only by the boundary conditions or 
any other relevant means of generating turbulence. However, he suggests that for high Reynolds 
numbers, when the turbulent motion is governed by a strange attractor, the symmetries of the 
governing equation are restored asymptotically for small scales. Velocity differences are used to 
localize the symmetry to small scales. 

The paradox inherent in this argument is that we cannot write governing equations for the 
velocity differences, exclusively in terms of velocity differences. A nonlinear term involving the 
velocity field, representing the sweeping interactions, is inevitable. As we shall argue below, the 
stronger homogeneity assumption u G 'Ki{A) is required to drop this term. Furthermore, we will 
argue that u G 3ii(A) is also required to derive the 4/5-law, which is the first step in Frisch's 
argument. Similar concerns were raised recently by Frisch [12] who questioned the self-consistency 
of local homogeneity, both in the sense of incremental homogeneity and in the sense of Kolmogorov. 

As for Kolmogorov, in his first paper [2], he defined local homogeneity in a very interesting 
way. Instead of using the Eulerian velocity differences w;a(x, x',t), he used the following quantity: 



Here, Y represents the approximate displacement of a fluid particle that is being used as a frame 
of reference. Because of its dependence on the velocity field, it is itself a stochastic variable. Kol- 
mogorov employed the probability density function of w^^; in his definitions. Furthermore, he 
included the requirement of local stationarity in his definition of local homogeneity. As will become 



Y(xo,to|x,t) = X - xo - (t - to)u(xo,to) 
Wi<roKxo,io|x,i) = u(Y(xo, io|x, i), - u(xo,to)- 



(24) 
(25) 
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apparent in section 4, Kolmogorov's representation of velocity differences is in fact a precursor of 
the quasi-Lagrangian representation, and we shall call it the Kolmogorov quasi- Lagrangian repre- 
sentation. Although Kolmogorov does not discuss explicitly the problem of sweeping interactions, 
it is interesting that he foresaw to this extent the need for an non-Eulerian representation of the 
velocity field. 

Another curious feature of the Kolmogorov definition is that it appears to use a conditional 
ensemble average conditioned on the statement u(xo,to) = v instead of the usual unconditional 
ensemble average, and includes independence with respect to v as part of his definition of local 
homogeneity. Equivalently, one may use a conditional average conditioned on the location of the fluid 
particle Y(xo, tQ\x, t) = y and assume independence with respect to y. The equivalence depends on 
using Kolmogorov's quasi-Lagrangian transformation, and is not applicable if one replaces it with 
the Belinicher-L'vov quasi-Lagrangian transformation. 

For the case t = to it is easy to see that the Kolmogorov definition is as strong as u G JCq. 
However, more generally, the velocity differences used by Kolmogorov are evaluated at two different 
times t and to- For this reason, I find it very unlikely that Kolmogorov's definition can be shown 
to be as strong as u G On the other hand, the use of the conditional ensemble average and the 
assumption that that average is independent of v probably strengthens the definition in unforeseen 
ways and may have some interesting consequences. In section 4.4 we show that a modified version 
of the Kolmogorov definition of local homogeneity is equivalent or stronger than 'Ki , thus providing 
the assumptions needed to prove the 4/5-law and eliminate the sweeping interactions. 

A detailed discussion and criticism of Kolmogorov's definition of local homogeneity is given by 
Frisch [12]. In the same paper, contrary to his prior work [7, 72], Frisch distinguishes the term local 
homogeneity from incremental homogeneity, and assigns Kolmogorov's definition as the definition 
of local homogeneity. We have seen that the definition of local homogeneity by Kolmogorov includes 
an assumption of incremental stationarity and also an assumption of a type of random Galilean 
invariance (i.e. independence with respect to v). Incremental stationarity can be true even when 
incremental homogeneity is not true. Furthermore, as I shall argue in section 5 of this paper, 
we should intend to derive random Galilean invariance from the theory instead of assuming it. 
Consequently, the original definitions [7, 72] of Frisch have the practical advantage of conveniently 
separating these assumptions from each other, and the conceptual advantage of not assuming too 
much. 



3.3 Balance equations and sweeping 



The clearest way to analyze the effect of the sweeping interactions on the theory of hydro- 
dynamic turbulence is by employing the balance equations of the Eulerian generalized structure 
functions. These balance equations were introduced by L'vov and Procaccia [50] in a landmark 
paper, and they arc derived as follows. 

The Navier-Stokes equations, where the pressure term has been eliminated, read 

+ T^pd^iUfsU^) = uV\a + '?a(sft3, (26) 

where CPq/j is the projection operator defined as 

= Sap - dadffV-^ (27) 

and da represents spatial differentiation with respect to x^. Repeated indices imply summation of 
components. The balance equations are obtained by differentiating the definition of F„ with respect 
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to time t and substituting the Navier-Stokes equations. This leads to exact equations of the form 



dt 



+ Dn = vJn + Qn, 



(28) 



where Dn represents the contributions from the nonlinear term, J„ the contributions of the dissi- 
pation term, and Qn the contribution from the forcing term. To write the terms concisely, we use 
the following abbreviations to represent aggregates of arguments: 

X=(x,x') 

{X}„ = {Xi,X2,...,X„} (29) 

{Xa}n = {Xl, . . . , Xfe_i, Xfc+i, . . . , X„}. 

The terms themselves read as follows. The forcing contribution is given by 



'?akli{fp{^k,t) - Jp{^k,t)) 



The dissipation term is given by 

n 

_^aia....a„(^X}„,i) = D„F„"i"--" ({X}„, i) = ^(V^^ + V^, JF„"1"--"«({X}„, t). 



(30) 



(31) 



k=l 



where V^j, differentiates with respect to xj., and similarly V^/^ differentiates with respect to x'y^. 

The remarkable result, shown in [50], is that the term D„ that represents the contribution of 
the nonlinear term can be rewritten as D„ = 0„F„+i + where 0„ is a linear integrodifferential 
operator with general form 



n+1 



^ ^ ^nk-^n 

k=l 

Gnk^n+l = J 0(Xfc, Yi, Y2) -F„+i({X}^, , Yi, Y2,t) dYitfY'2, 

and In is given by 



k=l \ 

where lX^({X}„,i) is defined as 
1 " 

U.a({X}„,t) = — ^ (na(Xfc,t) +na(x'fe,t)) . 



JJwai (X;,t) 



1=1 



(32) 
(33) 

(34) 
(35) 



k=l 



The first term, On^n+i, includes the effect of pressure and part of the advection term, and the 
detailed form of the operator 0^ has been given in Ref [76]. The second term, /„, represents 
exclusively the effect of the sweeping interactions. 

This decomposition makes rigorous the notion that the nonlinear interactions in the Navier- 
Stokes equations consist of local interactions that are responsible for the energy cascade and sweep- 
ing interactions which would disrupt the energy cascade if they contaminated the inertial range. 
It also exposes the conditions under which the sweeping interactions can be neglected. We learn 
that if the ensemble average of the velocity product that appears in the definition of /„ is invariant 
under a spatial shift, then the derivatives of that ensemble average will add up to zero. And here 
lies the heart of the problem. The assumption u G "Kq by itself is not sufficient to set In = 0. Global 
homogeneity u G ?{ is sufficient, but it is a stronger assumption than what is required. 
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Remark 2 The condition of incremental homogeneity, is written as u G %q{A). The homogeneity 
condition needed to eliminate the sweeping interactions over the domain A is u ^ 'Hi(A). 

It should be noted that the local term On-^n+i and the dissipation term J„ preserve the in- 
cremental homogeneity condition u G 5{o(-A)- The two terms in the balance equations that can 
potentially violate incremental homogeneity, are the sweeping term 7^ and the forcing term Qn- 
Asymptotic incremental homogeneity cannot be disrupted in the inertial range by the forcing term if 
the forcing spectrum is confined to large scales. The uncontrolled quantity is the sweeping term 
Recently, Prisch [12] questioned the consistency of local homogeneity, in the sense of Kolmogorov, 
and incremental homogeneity as a framework for studying hydrodynamic tTirbiilcnce. We see that 
incremental homogeneity can be a consistent framework on the condition that the sweeping term 
In is dominant only at large scales with its influence forgotten as the energy cascades to smaller 
scales. If that is the case, then none of the other terms in the balance equations violate incremental 
homogeneity. This is discussed in further detail in section 5. 

3.4 Remarks on the A/5-law proof 

In his second paper, Kolmogorov [3] employed an argument that is distinct from dimensional 
analysis to explain the claim that (^n = n/3. He derived the 4/5-law from which he obtained (s = 1, 
and used a scaling assumption to obtain (^2 = 2/3. Prisch's [7, 72] contribution was his observation 
that the scaling argument can be extended to account for all the scaling exponents C„. With this 
extension, Kolmogorov's second paper [3] is then an equivalent reformulation of the dimensional 
analysis argument of his first paper [2] . The superiority of the extended argument is that at least 
one of the scaling exponents is established rigorously. One also bypasses the universality criticism, 
appearently attributed to Landau, of the original similarity hypothesis of Kolmogorov. In his book, 
Frisch [7] gave a more detailed account of his argument, but he didn't derive the 4/5-law on the 
basis of incremental homogeneity and incremental isotropy as prescribed by his framework; he 
used instead global homogeneity and global isotropy. The same holds for the alternative proof by 
Rasmussen [77]. An old proof by Monin [78] and Monin and Yaglom [73] claimed to prove the 
4/5-law on the basis of local homogeneity and local isotropy, but it was criticized by Lindborg [79]. 
The criticism was addressed by Hill [80] who gave a corrected proof. 

In particular, the criticism of Lindborg [79] was that it was not proved that the correlations 
involving the pressure field gradient and the velocity field can be eliminated on the basis of local 
isotropy from the equation that governs the time derivative of the second order structure function 
tensor. Hill [80] resolved this objection by supplying the needed proof. The principle behind the 
proof is reflected, in a wider sense, by the mathematical form of the general sweeping term 
where there is only a local differential operator. The elimination of the nonlocal integral operator 
from /„ represents the elimination of any contributions by the pressure gradient term to /„ that 
would break incremental homogeneity. The pressure gradient does contribute to the term OnFn+i 
a non-local integrodifferential operator. However, because O^F^+i can be expressed exclusively in 
terms of the velocity differences, it preserves incremental homogeneity. 

Nevertheless, the proof by Monin and Yaglom [73], as far as our intentions are concerned, has 
an additional shortcoming, which has also been noticed independently by Frisch [12]: it concerns 
the elimination of the terms associated with the sweeping interactions. If we refer to the part of 
the discussion leading to equation (22.14) of Monin and Yaglom [73], we learn that they are using 
the Belinicher-L 'vov quasi- Lagrangian transformation to eliminate the sweeping interaction terml 
This can be made more clear if the reader compares the argument involving the two unnumbered 
equations that precede equation (22.14) of Monin and Yaglom [73] with section 4 and appendix A. 
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The intention of this argument, according to Monin and Yaglom [73], is to "... transform the Navier- 
Stokes equations so that they contain only the velocity differences and their derivatives". This is 
precisely what the quasi-Lagrangian formulation does. As we shall argue in the next section of this 
paper, applying the inverse transformation back to the Eulerian representation uses an assumption 
of homogeneity stronger than incremental homogeneity, but this time u G IK* . 

As far as the theory of the scaling exponents is concerned, it is only necessary to know ^3. An 
elegant way to calculate (^3 is from the solvability condition of the homogeneous equation O2-F3 = 
[48, 56]. The idea here is to use the conservation of energy to show that 

/ ^ /X ^ djSajru) - S3{ri2')] , 1 d[S3{rr2') - S3{ri>2)] 

U2-r3lXl,X i,X2,X 2j = h 



2 dri 2 dry (36) 

C3- 

1'2 



_ 4r„C3-i _ C3-1 I „C3-i _ „C3-i 
~ ^l' 12 '12' ^ ' 1'2' ' 1' 



where ri2 = ||xi — X2II, etc. It follows that the equation O2-F3 = will be satisfied for any config- 
uration of velocity differences if and only if ^3 = 1- The homogeneous equation can be obtained 
from the balance equations in the limit of infinite Reynolds number. For the case of finite Reynolds 
number, there is a homogeneous and particular solution to the generalized structure functions that 
are linearly superimposed [81, 82]. Then the calculation of ("3 is relevant only for the homogeneous 
solution. Aside from this issue, this argument too requires that we set I2 = 0. Dropping I2 cannot 
be justified under incremental homogeneity, in the sense of u G 5{o) and it requires the condition 
u € Jfi(yi). We arrive then to the following conclusion. 

Remark 3 The homogeneity condition needed to establish Cs = 1 over the domain A is u & "Ki (A) . 

It should be noted that even though Hill [80] has claimed to show the 4/5-law on the basis of 
local homogeneity and local isotropy, his definition of local homogeneity is mathematically stronger 
than the definition u £ Jfo(^) used in the Frisch framework, and it is in fact very similar to 
u G 'Ki{A) (also see section 4.1 of [83]). Consequently, while his proof correctly follows from his 
stated assumptions, it cannot be used from within the Prisch framework of hypotheses to prove the 
4/5-law without invoking additional assumptions. 

It is possible to derive a rigorous version of the 4/5-law that does not require assumptions 
of homogeneity, isotropy, stationarity, and not even an ensemble average [84-86] . This is done by 
rephrasing the statement to be proven. Specifically, it has been shown that 

, , /•*+^* , f dn(A) f dx S3(x,rAe) 4 
lini lini lim / dr / — ^ / ' ^ = --£3, 37 

M^Or^Ou^oJt Jso(3) ^TT ^3 ^(S) r 5 

for almost every (Lebesgue) point t in time, where e is a unit vector, S C is a local region in 
a periodic boundary domain (topologically equivalent to a torus) with volume V{'B), and £3 is 
the local dissipation rate over the region 23 given by 



£3 = lim 



^y"dx£(x,0, (38) 



where e(x, t) = (l/2)z^ {sa,/3{x,t)sap{x,t)) is the dissipation rate density at (x, t) and s^fs = dctUjj + 
djjUa is the local strain tensor. A similar result was obtained earlier by Nie and Tanveer [87]. 

It should be noted that this result does not contradict our previous remark. Although the need 
to make assumptions appears to have been eliminated, this is done so at the price of proving a 
statement that is mathematically weaker. In the original formulation of the 4/5-law, aside from 
an ensemble average, all the integrals are absent. These integrals represent an interesting way of 
obviating the symmetry assumptions needed to prove the 4/5-law in its original formulation. 
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Recently, there has been considerable interest in extending the 4/5-law to account for deviations 
from the theoretical prediction caused by the violation of incremental isotropy [88-92]. From the 
viewpoint of the experimentalist these extensions make it possible to confirm the validity of the 
4/5-law against experimental data. Prom the viewpoint of the theorist, deviations from incremental 
isotropy can be accounted for with the SO{3) group decomposition method [56, 57]. 



4 The qucisi-Lagrangian formulation 

The essence of the quasi-Lagrangian formulation (also called the Belinicher-L'vov transfor- 
mation) is to look at turbulence using a fluid particle as a non-inertial frame of reference. The 
representation is Lagrangian because we involve fluid particles, but it is not completely Lagrangian 
because the fluid particle trajectory is only used to define a new frame of reference, and we continue 
to look at the velocity field in an Eulerian manner. It is understood, of course, that the only in- 
teresting statistics are those involving points within a sphere centered on the moving fluid particle 
with radius on the order of the integral length scale ^q. 

Let Mq(x, t) be the Eulerian velocity field, and let pa{xQ, to|i) be the position of the unique fluid 
particle initiated at (xo,io) at time t relative to its initial position at time tQ. The transformation 
is done in two steps. First, we introduce Va{^o,to\:x.,t) as the Eulerian velocity with respect to the 
original inertial frame of reference with a space shift that follows the fluid particle: 



J to 



(39) 

Wa(xo,to|x,t) = Ua(x + /9(xo,to|t),t). 

Then, to complete the transformation we must subtract the velocity of the fluid particle uniformly, 
so that the particle itself will appear to be motionless: 

d 

W^a(xo,io|x,i) = Va(xo,io|x,t) - ^Pa(xo,to|*) = Va{^,to\x,t) - Va{^,to\:s.o,t) ^^^^ 
= Wa(x + /9(xo, to [*)>*) -^a(xo + p(xo,to|i)>*)- 

We deflne Wai'x.o, *o|x, t) as the quasi-Lagrangian velocity field, and introduce the quasi-Lagrangian 
velocity difference Wq;(xo, to|x, x', given by 

WQ,(xo,to|x,x',t) = Wa{xo,to\x,t) - ttJQ, (xq, to |x', t) = f q, (xq, to |x, t) - f q, (xq, to |x', t) . (41) 

Differentiating with respect to time, and substituting the Navier-Stokes equations, gives an equation 
of the form 

^ + Vaf^jWpW^ = + Vl,)Wa + Fa , (42) 

where ^^(xo, to|x, x', t) is the quasi-lagrangian forcing, and VaP'y is a bilinear integrodifferential 
operator of the form 

Va^-yWpW^ = JJ dXpdX^ VaPji^o\y^a, X^, X^)Wp{Xp)W^{X^), (43) 

with y(xo|XQ, X^, X-y) the corresponding kernel (see appendix A for more details). The remarkable 
feature of this equation is that all the terms, and most especially the nonlinear term, are written 
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in terms of velocity differences. Fundamentally, this is the reason why the quasi-Lagrangian trans- 
formation eliminates the sweeping interactions and renormalizes the MSR diagrammatic theory. 

The key issue is whether it is possible to switch back to the Eulerian representation without 
reintroducing the sweeping interactions. In a short appendix to their paper, L'vov and Procaccia 
[13] showed that in stationary turbulence the ensemble average of the same time quasi-Lagrangian 
velocity differences is equal to the ensemble average of the corresponding Eulerian velocity differ- 
ences. The same appendix is also found in a previous unpublished paper [93]. The proof requires 
stationarity of the Eulerian velocity field, and incompressibility. A homogeneity condition is also 
used, which is described as "translational invariance" . 

In this section, we would like to carefully re-examine this proof, the assumptions needed to make 
it work, and the relationship between this result and other claims that one might reasonably make 
about the quasi-Lagrangian velocity differences. Part of our motivation is the crucial importance of 
this result; the entire L'vov-Procaccia theory [13, 46, 47, 51-55] stands or falls on the validity of this 
argument. As we have discussed previously in section 2, the L'vov-Procaccia theory and the Giles 
theory [64] are the only two theories that can explain mathematically the reason why the inertial 
range of three-dimensional turbulence has intermittency corrections. Our main interest is to show 
that the proof requires that we assume u G Jf* , which is a stronger condition than what is actually 
needed to eliminate the sweeping interactions or to prove the 4/5-law (u G Jii{A)). Preliminaries 
are given in section 4.1 and section 4.2, and the proof itself is discussed in section 4.3. 



4-1 Characterizations of the claim 



Let 9'„(xo, to|{x, x'}„, t) be the generalized structure function in the quasi-Lagrangian repre- 
sentation, defined as 



9'n(xo,to|{x,x'}„,t) 



.k=l 



(44) 



The claim of L'vov and Procaccia [13] was that it can be shown that 



^„(xo,to|{x,x'}„,t) = F„({x,x'}„,t),Vn G r 



(45) 



which can be rewritten equivalently as 

Wa(xo,to|x,x',t) U^a(x,x',t). (46) 



As a first step, consider the following easy-to-prove propositions which give equivalent characteri- 
zations of the claim (46): 

Proposition 1 The claim (46) holds if and only if the quasi- Langrangian velocity is incrementally 
stationary with respect to to: 

Wa(xo, to + At|x, x', t) Wa(xo, to|x, x', t), VAt G M - {0} (47) 
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Proof: (=^): Assume that the claim (46) holds. Then, it follows that 



Wa{xo,to + At\ji,x',t) ~ u;a(x,x',i) ^ Wa(xo,io|x,x',i). (48) 

{<=): Now assume that the quasi-Lagrangian velocity field is incrementally stationary with respect 
to to- Using the evaluation lim.tg-^t Pa(xo,io|*) = 0, it follows that 

/ X / X x' / 

Wa(xo,to|x,x ,t) ~ lim VFa(xo,to|x,x ,t) ~ lim u;a({x, x } + p(xo, ioK)> (49) 
to— >t to—*t 

^'u;«(x,x',t), (50) 



X5X 



□ 

Proposition 2 Assume incremental stationarity on the Eulerian velocity field: 

Waix,x',t) Wa(x,x',t + At), VAt G R. (51) 

Then, the claim (46) /lo/ds if and only if the quasi-Langrangian velocity is incrementally stationary 
with respect to t: 

Wa(xo,to|x,x',t) ~ VF«(xo,to|x,x',t + At), VAt G M - {0} (52) 
Proof: (^): Assume that the claim (46) holds. Then, 

Wo,(xo,to|x,x',t + At), i(;c,(x,x',t + At) -u;q(x, x', t) VFq(xo, to|x, x', t). 

(-^): Now assume that the quasi-Lagrangian velocity field is incrementally stationary with respect 
to to. Using the evaluation limt(j_>t Pq,(xo, to|t) = 0, it follows that, 

WQ,(xo,to|x,x',t) lim VFQ,(xo,to|x,x',t) lim Wa({x,x'} + p(xo, to|t), t) 
t— >to t— >to 

= u;a(x,x',to) '''^ u;a(x,x',t). 

□ 

The implication of these propositions is that the relationship between the Eulerian and the 
quasi-Lagrangian formulations can be established as an immediate consequence of stationarity of 
the quasi-Lagrangian velocity field with respect to to- As Lebedev and L'vov [93] noted, the variable 
to does not appear anywhere in the quasi-Lagrangian Navier-Stokes equations (42), consequently the 
form of the governing equations allows stationary solutions with respect to to- However, to assert 
that the quasi-Lagrangian velocity field is stationary, it is necessary to assume that the quasi- 
Lagrangian forcing field is also stationary. Since the definition of the quasi-Lagrangian forcing field 
entangles the Eulerian forcing field fa with the trajectory field pa, and since the trajectory field 
itself is not time invariant (due to the initial condition Pq,(xo, to|to) = 0), we cannot make this 
assumption without justification. This was the reason, cited by Lebedev and L'vov [93], for the 
rigorous proof which is the topic of this section. Proposition 2 shows that assuming stationarity in 
the quasi-Lagrangian representation is sufficient to prove the claim (46) , and thus this assumption 
implicitly introduces in the Eulerian frame the conditions needed to prove the claim. 
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4-2 MSR theory for Lagrangian trajectories 



The governing equation for pa is 
dpa{yio,to\t) 



at 



Ua{xo + pixo,to\t),t), (53) 



with initial condition pa{xQ,to\to) = 0. Deriving the stationarity condition (47) requires an MSR 
theory where the velocity field Ua can be thought of as the forcing field with known statistical 
properties, and the Lagrangian trajectories field pa as the governed field whose properties wc wish 
to deduce. Unfortunately we may not apply the standard MSR theory because the equation itself 
does not assume the standard form Na[p\ = Ua with Ua independent of p and furthermore the 
initial condition is set at a finite time to and not at to —oo. We need to develop the statistical 
theory from scratch, and for that purpose the path integral formulation is most expedient. 

Note that every value of to corresponds to a distinct initial value problem. We may therefore 
treat, for the purposes of the statistical theory, the field pa as a function only of xo,t and let to 
to be taken at a fixed value. We can also go a step further and note that for every value of xo 
the governing equation is an ordinary differential equation. It follows that in constructing an MSR 
theory for pa we have two options: Wc may construct a statistical theory for the restricted problem 
in which xo is also fixed and the field pa, is taken as a function of only t, or a theory for the full 
problem in which only to remains fixed and pa is taken as a function of xq and t. In the restricted 
case we cannot calculate correlations between fields pa with different values of xo. In the full case, 
we can. For our needs, the restricted statistical theory will be sufficient. 

Introduce an operator Qxg [p] via the kernel 

Qa^[p](*,y,^) = ^ap5{t - T)5{y - xo - p(t)), (54) 
such that 

//•+00 
dy I dr Q'^°Jp\{t,y,T)up{y,T) = Ua{yio + p{t),t). (55) 
J-oo 

We also introduce a functional 'C/xo,toM ^^^^ constructs Pa from the velocity field. This operation 
is of course admissible in both the restricted and the full theory. 

Since p = X^xo.toM equivalent to pa = Qxo[p]^ia, it follows that there exists a functional J[u] 
such that 6[pa — Qxo[p]ua\ = J[u]S[p — ■^xojtoM] which can be evaluated by integrating both sides 
over p: 

J[u] = [ Vp6[pa-Q^o[Pha]- (56) 

The integral is a Feynman path integral [94] (a pedagogical introduction is given in Ref. [95]). Here, 
J'(to) is the domain of integration and it is defined as the set of all Pa(t) that satisfy the initial 
condition Pa(to) = 0. Wc also define T as the set of all possible paths. Suppose we would like to 
evaluate the ensemble average {M[u,p]) where M is some arbitrary functional of pa and Ua- We 
treat the velocity field Ua as a forcing field with known statistics. We assume then that we know 
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how to evaluate the ensemble average of any expression in terms of the velocity field. We have: 

(M^p]) = {M[u,^:^,,toH) = ( [ 'D0M[u,g]d[g-i:^„toM) 

\Jnto) I 

Vg {M[u,g]5[g-ii^,,to[y]]) 



L 



! 'Dg {M[u, g]J-\u]5[ga - Q^Ua]) 



[ / 2)/3 expiiPaQa){M[u,g]J-^[u]exp{-if3aQ^o[Q]ua))- 

J nto) Jv 



into) 

Here, wc have used the convention that repeated Greek indices imply integrating temporal coordi- 
nates throughout their domain in addition to summation of vector components. For example, the 
expression Pa Pa is an implicit abbreviation for 

/+0O 
dt pa{t)Pa{t). (57) 
-oo 

We also use the formal representation for the delta functional 

6[u] = [ Df3 expiiPaUa), (58) 

which is valid in the sense of generalized functional distributions. 

If the velocity field is incompressible, it can be shown that J[u] = 1. A detailed proof of this 
result is given in appendix B. Then, the stochastic theory simplifies to: 

{M[u,p])= [[ DgDp exp{iPaga){M[u,g]exj>{-iPaQ^,[g]ua)). (59) 

J Jnto)xv 

This statement is a concise expression of the statistical theory for Lagrangian trajectories. 
4-3 Transform back to Eulerian representation 

We now use the statistical theory to derive the relationship between the quasi-Lagrangian 
correlation and the Eulerian correlation. The proof given here follows the one given by L'vov and 
Procaccia [13], but it is presented in more detail to show the underlying assumptions. The proof 
makes an essential use of Eq. (59) derived above. 

The argument is essentially based on the following identity: 

Pairo,to + At\t) = pairo,to\t) - pa{ro,to\to + At). (60) 

To see why this is true, note that the expression on the right-hand side satisfies the governing 
equation for the + At problem, and it also satisfies its initial condition. Therefore, by uniqueness, 
the right-hand side has to be equal to the left-hand side. 
To facilitate with calculations, we define Mxq [P, p] as 

Mxo [P, p] = {M[u, p] exp(-i/3«Qxo [pM) (61) 
V[u,p]exp (^-i J dtpa{t)ua{^o + p{t),t)^y (62) 
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and we also use the notation M(xo,to) for the ensemble average {M[u,p]) evaluated under a given 
choice of xq and to. Consequently, we may write 

M(xo,io)= // DpD/3exp(i/3„/9„)Mxo[/3,p] (63) 
= [[ 2)pD/3exp (i I dt I3^{t)^^\ Mx„[/3,p]. (64) 

J J y (to ) X y \ *f / 

The key statement to be proven is the following proposition, that shows the connection between 
stationarity in the quasi-Lagrangian representation and homogeneity in the Eulerian representation. 
Proposition 3 // the. velocity field u is incompressible, then 

V/3 e T : Vx G M"' : Mxq [/?, P + x] = Mxo [P, p] =^ VAi G M : M(xo, to + At) = M(xo, to) (65) 



Proof: 

To facilitate our argument, introduce a new field Aq, defined as equal to the right hand side of 
(60). 

Xa{ro,to\t) = Pa(.ro,to + At\t) (66) 

= Pa{ro,to\t) - Pa{ro,to\to + At) = Ba/3P/3 ■ (67) 

The connection between Aq, and pa is linear, in the sense that we can construct an appropriate kernel 
Bai3 made of delta functions that transforms one field into the other. The functional determinant 
of B is equal to 1, so a change in variables under the path integral does not introduce an additional 
factor, namely DA = Dp. This is usually true with simple transformations, such as space shifting and 
rotations, because they merely reshuffle the order in which we integrate over all possible histories. 
In this case, we need to take into account that the permissible histories are constrained by the 
initial condition Aa(ro, tol^o + ^t) = which is different from the initial condition of the field pa- It 
follows that, while there is no need to introduce a functional determinant, the domain of integration 
has to change from J'(to) to y(to + At) . Finally, it is easy to see that dXa/dt = dpa/dt and the 
hypothesis implies that Mxo[/?, p] = ^xoiP, A]. We may then write: 

M(xo, to)= [[ DpD(5 exp (i [ dt p^{t)^^) Mxo p] (68) 
JJv(to)xV \ J at ) 

= 11 DAD^exp [i [ dt(3a{t)^^) M^,[/3,X] (69) 

J J 7(10+ At) \ J at J 

= M(xo , to + At) VAt G R. (70) 

□ 

Proposition 4 If Ua is incompressible, then 

u G 5£* ^ VAt G M : Wa{^o,to + At|x,x',t) T^a(xo,to|x,x',t). (71) 
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Proof: Let n G N* be given, and define the functional M[u,p] as 



M[u,p] = Yl W"afc(xo,to|Xfc,x'fc,t) 

k=l 
n 

k=l 

Consequently, the functional Mxq p] reads 



(72) 
(73) 



Mxo [P, P] = (^M[u, p] exp i^-i J dt Pa{t)Ua{yiO + p{t), t) 



(74) 



m=0 
X 



n^A(^o,io|xi,x'i,t) 



.1=1 



JJ (xo + p(tjfc),ifc) 



.fc=i 



(75) 



From the assumption u G IK* we see that the ensemble average in the equation above is invariant 
with respect to a uniform spatial shift. It follows that Mx,, [/?, P + x] = Mxo[/3, p],Vx € , and 
using proposition 3, this implies that M(xo, to + ^t) = M(xo, ^o); VAt G M. Consequently, we have 



^•^(xctoKXIn,*) =M(xo,to) =M(xo,to + At) 

= :J„(xo,to + At|{X}„,t), Vn e N*. 



(76) 
(77) 



□ 

The claim (46) follows by combining proposition 4 with proposition 1. It should be noted 
that once the relationship between quasi-Lagrangian correlation functions and Eulerian correlation 
functions is established, it can be easily extended to response functions as well without making 
any further assumptions. Starting from the stationarity condition (51), we deduce from the quasi- 
Lagrangian formulation of the Navier-Stokes equations that the quasi-Lagrangian forcing field is 
also stationary. Then, we may use an MSR theory on the quasi-Lagrangian Navier-Stokes equations 
to obtain stationarity on the response functions. From there, the relationship between the quasi- 
Lagrangian response functions and the Eulerian response functions can be easily established. 



4-4 A derivation via conditional local homogeneity 

The artifact introduced by the quasi-Lagrangian formulation is that the turbulent velocity 
field is being perceived from the viewpoint of an arbitrary fluid particle whose own motion is also 
stochastic. Consequently, to relate the quasi-Lagrangian correlation tensor 9'„(xo, to|{X}ri, t) with 
the Eulerian correlation tensor i^„({X}„, t) a certain sense of homogeneity is required to ensure that 
the velocity field is being perceived by the fiuid particle in the same way regardless of the actual 
position of the particle. Our analysis of the proof, given previously, has shown that the homogeneity 
condition used by the proof is stronger than the condition u G JCi required to eliminate the sweeping 
interactions. What is particularly interesting about the stronger condition u G IK* is that it requires 
translational invariance from a group of the many-time correlation tensors -F^„. 
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Let us now consider an alternative approach. Introduce the conditional correlation tensor de- 
fined as 



5'n(xo,io,y|{X}n,t) = ( Jl T4^a;^(xo,to|Xfc,x'fc,t) 

\a:=1 
/ n 



pixo,to\t)=yj (78) 
p(xo,io|t) =yV (79) 



\k=l 



This definition is identical to the definition of the quasi-Lagrangian correlation tensor 5'n(xo, to|{X}„, t), 
except that the ensemble average is replaced with the conditional average predicated on the fluid 
particle being located at position y at a given time t. Let p(xo, to\x, t) be the probability that a fluid 
particle originating at (xq, to) will be located at x at time t. It follows that the Eulerian correlation 
tensor F„({X}„,t) and the quasi-Lagrangian tensor ?'„(xo, to|{X}„, t) are given by 

F„({X}„,0 = ydy:J„(xo,to,y|{X}„-y,Op(xo,io|y,t) (80) 

^„(xo,to|{X}n,t) = y"dy^„(xo,to,y|{X}„,t)p(xo,to|y,i)- (81) 

It is trivial to see that if ?'„(xo, y|{X}n, t) is invariant with respect to {X}„ t— »■ {X}„ + Ax 
(conditional local homogeneity), then the Eulerian correlator and the quasi-Lagrangian correlator 
will be equal. 

Now let us consider the implications from a Kolmogorov-like definition of local homogeneity 
where wc assume that J'„(xo, to) y|{X}„, t) is independent of y without assuming invariance with 
respect to {X}„ i-^ {X}n + Ax. Then we have 

^„(xo,to|{X}„,t) = y"dy^„(xo,to,y|{X}„,t)p(xo,to|y,t) (82) 

= 5'„(xo,to,y|{X}„,t) J dy pixo,to\y,t) (83) 

= 3^„(xo,io,y|{X}„,i), (84) 

and 

Fn{{X}n,t) = y"dy:J„(xo,io,y|{X}n-y,t)p(xo,to|y,i) (85) 

= J dy ?'„(xo,to|{X}n -y,t)p(xo,io|y,i)- (86) 

To establish equality between the Eulerian correlation tensor F„({X}„, t) and the quasi-Lagrangian 

tensor 9"„(xo, to|{X}„, t), we also need 3"„(xo, to|{X}n, i) to be invariant with respect to {X}„ 
{X}„ + Ax. Unfortunately, this can only be established if we also assume an ergodic-likc hypothesis 
that p(xo, to|x, t) is independent of x. It is reasonable to expect this hypothesis to hold for t 3> to. 
Then, it follows that 

3^„(xo,io|{X}„ + Ax,t) = ydy3^„(xo,io,y|{X}„ + Ax,i)p(xo,io|y,i) (87) 

= J dy 9'n(xo, to, y + Ax|{X}„ + Ax, t)p(xo, to\y + Ax, t) (88) 
= J <^y 3'n(xo,io,y|{X}„,i)p(xo,io|y,*) = 3'n(xo,io|{X}n,t). (89) 
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This result, combined with Eq. (86) and (84), imphes that 



9^„(xo,to|{X}„,t) = i^„({X}„,t) 



(90) 



which in turn, combined with Eq. (89) also gives the u G !Ko condition: 
F„({X}„ + Ax,t) = F„({X}„,t), VAX G R". 



(91) 



It appears that the ergodic assumption on p(xo,to|x,i) is necessary to derive Eqs. (89), (90), and 
(91). The assumption of invariance of 3"„(xo, y|{X}n, i) with respect to y is not sufficient. 

As we have mentioned previously, in his first paper, Kolmogorov [2] also defined local homo- 
geneity using a conditional ensemble average conditioned on the fluid velocity at the reference point. 
Because of the approximate nature of the quasi-Lagrangian transformation used by Kolmogorov 
(which is not identical to the quasi-Lagrangian transformation of Belinicher and L'vov), his def- 
inition can be rephrased in terms of a conditional average on the location of the fluid particle. 
We have shown that if one uses the Belinicher-L'vov quasi-Lagrangian transformation instead of 
the Kolmogorov quasi-Lagrangian transformation, and changes the conditional ensemble average 
from using the velocity of the reference fluid particle to using the location of the reference fluid 
particle, then this modified definition of local homogeneity combined with a reasonable ergodic-like 
hypothesis does eliminate the sweeping interactions. We may conjecture that Kolmogorov had the 
elimination of sweeping in mind when he formulated his definition, but there is no such explicit 
indication in his papers. 

5 How the elimination of the sweeping interactions should be justified 

We have seen that when using the quasi-Lagrangian transformation we end up making the 
homogeneity assumption u G IK^ which is much stronger than the assumption of incremental 
homogeneity u G ^o{A) of the Frisch framework of hypotheses, otherwise we cannot return back 
to the Eulerian representation. Furthermore, this homogeneity assumption is introduced implicitly 
just by assuming stationarity in the quasi-Lagrangian representation, even if we don't wish to go 
back to the Eulerian representation (see proposition 2). The question that we would like to consider 
now is whether the utility of the theoretical work [13, 46, 47, 51-55] that relies on the transformation 
itself is diminished. We would like to claim that this is not the case, and define a line of investigation 
that can clarify this further. From an experimental standpoint, the very existence of a robust energy 
cascade indicates that the sweeping effect is confined to the large scales, and therefore it can be 
neglected with impunity. The main question that needs to be addressed is: How should one justify 
theoretically the elimination of the sweeping interactions? 

5. 1 Elimination of the sweeping interactions 

It is widely accepted that the behavior of the structure functions in the inertial range does 
not depend on the statistical properties of forcing, as long as the spectrum of the forcing term is 
confined to large length scales. In a sense, as the energy cascades toward smaller length scales, the 
characteristic features of the forcing term are "forgotten" . One may conjecture that the sweeping 
interactions behave in a similar way as a large-scale forcing term whose effect is forgotten in the 
inertial range. We may base this conjecture on the fact that even though the required homogeneity 
symmetry u G 3ii(yi) may not hold exactly, it can be expected to hold asymptotically at small 
scales. Consequently, even though we cannot set /„({X}„) = exactly, we might expect this term 
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to b econie rapidly snicill wliGii tlic £iv6rcLg6 SBpciTcition R between the points {Xj-^^ goes to zero. 
But does it vanish rapidly enough? A rigorous argument would have to estimate how fast /„, as a 
function of R, is approaching zero in the small-scale limit R/io — > 0, and then calculate the scaling 
exponent A„ associated with the ratio 

{OnkFn+immn) J ' ^ ' 

where R is the scaling parameter and Iq the forcing scale. Then, provided that one starts with 
the assumption u € proving A„ > is also a proof that u € 'Ki{A) which is sufficient to 
eliminate the sweeping interactions. Let A„ be the scaling exponent of /n(-R{X}„). If we assume 
that the generalized structure functions F„(i2{X}„) satisfy the fusion rules [48, 50], then the scaling 
exponent of 0„fcF„-|-i(i?{X}„) is Cn+i ~ 1 '^^d it follows that A„ = A„ — (Cn+i — !)• The challenge, 
then, is to calculate the scaling exponents An which are not likely to be universal. 

It is easy to see that this argument cannot be extended to the inverse energy cascade of two- 
dimensional turbulence. In that case, the forcing term operates at large wavenumbers. Given that 
we can reasonably assume that the inverse energy cascade is local, we expect that the forcing 
term is forgotten in the inertial range. The problem is that the energy is now going towards 
small wavenumbers. As we have noted, J„ essentially measures how much homogeneity is violated 
at a given length scale R. At large length scales, the flow will begin to sense the violation of 
homogeneity caused by the boundary conditions which will in turn make the sweeping term I^i 
larger in magnitude. If it becomes comparable to the terms OnkPn+i, it will probably disrupt the 
inverse energy cascade. 

Numerical simulations have shown that it is possible to obtain an inverse energy cascade under 
certain conditions [96], but it can also be disrupted under other conditions [97-100]. Physically, 
this disruption arises from the spontaneous generation of long-lived coherent vortices that carry a 
significant amount of enstrophy. An explanation of this effect was given by Boffetta et al [96] , in 
terms of the "bottleneck" effect [101]. The general idea is that the behavior of the energy spectrum 
in the inertial range is modified at wavenumbers near the dissipation range because some of the triad 
interactions at these length scales are disrupted by the dissipation term, thus making the transfer 
of energy less efficient. It is reasonable to anticipate the same effect in a high-resolution simulation 
of the inverse energy cascade, where the cascade has manifested successfully, without being arrested 
by coherent structures. However, we would like to suggest that the deviations observed by Danilov 
and Gurarie [98-100] and Borue [97] are more likely to be caused by a similar effect where the triad 
interactions are disrupted hy the sweeping term In rather than the dissipation term DnFn at large 
scales. 

The coherent structures that appear in two-dimensional turbulence can be conceptualized as 
concentrated small blobs of very high vorticity that raise a two-dimensional "hurricane" in the 
velocity field around them. Thus, from an intuitive standpoint, it is reasonable to expect that their 
presence in the flow should amplify the sweeping effect. Prom a theoretical perspective, one can say 
that /„ excites a "particular" solution of the statistical theory for the correlators F„ which combines 
linearly [102, 103] with the "homogeneous" solution of the homogeneous theory OnFn+i = that 
corresponds to the inverse energy cascade. It follows then that to obtain an inverse energy cascade 
in the forced-dissipative setting, one requires a dissipation term at large scales which will not 
only dispose of the incoming energy, but will also damp out the sweeping term /„ over the entire 
range of length scales where it is comparable to OnkFn+i- It also follows that there should be 
a conspicuous discrepancy between the energy spectrum in the quasi-Lagrangian representation 
and the energy spectrum in the Eulerian representation, when the coherent structures provide 
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the dominant contribution to the Eulcrian energy spectrum. It is well-known that in the Eulerian 
representation the coherent structures contribute a dominant A;~^ term [97]. On the other hand, in 
the quasi-Lagrangian representation, one should recover the k~^^^ contribution from the underlying 
inverse energy cascade. It has already been established that the underlying k~^/^ spectrum can be 
recovered if the coherent structures are artificially removed, either by a wavelet technique [104], or 
more crudely [97,99]. If our conjecture holds, then it should be possible to obtain the same effect 
simply by transforming into the quasi-Lagrangian representation. 

We have referenced the inverse energy cascade of two-dimensional turbulence as an example 
where it is not safe to "eliminate" the sweeping interactions. The criticism of the quasi-Lagrangian 
formulation by Mou and Weichman [105] is essentially that it has not been demonstrated that it is 
"safe", in the same sense, to eliminate the sweeping interactions in the downscale energy cascade 
of three-dimensional turbulence. 

5.2 Alternatives to Lagrangian transformations 

It is possible to use the theoretical work based on the quasi-Lagrangian transformation in a way 
that requires only the assumption u € J£i (^l) instead of u G CK* . This can be done via the following 
line of argument: The quasi-Lagrangian formulation modifies the Navier-Stokes equations by re- 
defining the material derivative (see appendix A). The modified equation remains mathematically 
equivalent to the Navier-Stokes equation because the velocity field is reinterpreted from an Eule- 
rian field into a quasi-Lagrangian field. It is precisely this reinterpretation which necessitates the 
stronger assumption u G IK* to enable a return back to the Eulerian representation. On the other 
hand, if we accept the hypothesis that the sweeping interactions can be absorbed into the statistical 
forcing term, we can modify the equation of motion in precisely the same way without changing the 
interpretation of the velocity field. From there, one can derive the same balance equations (28) ex- 
cept that one will have In = 0, and consequently the only assumption that is being made implicitly 
is just u G IKi(yi). One may then proceed from this Eulerian modified Navier-Stokes equation and 
develop the L'vov-Procaccia theory [13, 46, 47, 51-55] with impunity, since the modified governing 
equation would have the same mathematical form as the quasi-Lagrangian Navier-Stokes equation. 

In geometrical language, the difi^erence between the quasi-Lagrangian transformation and what 
I propose is the following: In the quasi-Lagrangian transformation we perceive the flow from the 
point of view of a single fluid particle. The claim to be established is that the one-time statistical 
properties of velocity differences should remain invariant when switching between the inertial frame 
of reference and the non-inertial frame of reference defined by the fluid particle. My suggestion is to 
consider a transformation where for a given point in space and time (x, t) we perceive the flow from 
the point of view of whatever fluid particle just happens to be there at (xo,t). This leads to the 
Navier-Stokes equations for the Eulerian velocity differences Wq,(x, x',t). Then, the homogeneity 
assumption u € 3ii is sufficient to establish the claim that the one-time statistical properties of 
velocity differences Wa(x, x', t) will remain invariant under a transformation from the inertial frame 
of reference to a non-inertial frame of reference defined by the fluid particle at (xo,t). This claim 
is in fact mathematically equivalent to the condition /„ = which follows from u G !Ki. 

In connection with this argument, it is interesting to note that the idea of just modifying the 
Navier-Stokes equation was considered by Kraichnan [25] in 1964, who suggested a more crude 
modiflcation. This modification brute-forces locality in Fourier space by discarding triad interac- 
tions across a wide wavenumber interval and retaining only the local triad interactions. From the 
same paper we learn that Kraichnan suspected that there was a relationship between the quasi- 
Lagrangian transformation of Kolmogorov and the general idea of modifying the Navier-Stokes 
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equation in such a way but noted that bringing that out rigorously is difficult. In my view the 
quasi-Lagrangian transformation of Belinicher and L'vov, which is different from the Lagrangian 
transformation used by Kraichnan in his theories, is the key to finding possibly the best way to 
modify the Navier-Stokes equations in the way that Kraichnan and Kolmogorov intended. 

An alternative argument that was proposed by Yakhot [43] and used by Giles [64] to calculate 
a perturbation expansion for the scaling exponents Qn eliminates the sweeping interactions by 
modifying the statistical theory itself. This is different from the quasi-Lagrangian formulation and 
my proposal where the change is made on the governing equation and then propagated into the 
statistical theory. Again, to justify why one can modify the statistical theory requires the assumption 
u G JCi [A) or an argument justifying the hypothesis that the sweeping interactions can be modeled 
as large-scale stochastic forcing, which brings us back to the challenge of showing that A„ > 0. 



5.3 Estimating the scaling exponent A, 



The problem of calculating the scaling exponents A„ and A„ needs to be investigated primarily 
with numerical simulations and the analysis of experimental data. However, it is possible to make a 
speculative theoretical calculation, if we are willing to commit the following crimes against reality: 
First, we assume that the velocity field Uq,(x, t) can be modeled as a random gaussian delta- 
correlated (in time) stochastic field acting at large scales. Furthermore, we assume that the velocity 
field ?Xa(x, i) has an effect on the velocity differences ^^(x, x',t) via the sweeping interactions, but 
completely disregard the reverse effect of the velocity differences on the velocity field via eddy 
viscosity, and the fact that Uaiy^^t) and Wa{x,x.' ,t) are obviously constrained by the definition of 
Wai^jx' ,t). In other words, we assume that ^^(x, x',t) is advectcd as a passive scalar by Uq(x, t) 
and that tta(x, t) can be assumed to be a random gaussian delta-correlated in time field. Note that 
WcK(x,x',t) is still also forced by fa- 

We have shown in Appendix C that under these assumptions the sweeping term 7„({X}„,t) 
can be decomposed into three contributions: 

In{mn,t) = 4,(i)({X}„,t) +7„,(2)({X}„,t) +7„,(3)({X}„,t), (93) 

which are given by 

7;^(^--"({X}„,t) = {Up{{XU,t)) H-^----r^^{{XU,t), (94) 

n n 

i::sr''imn,t) = Y.I1 ^^:-r'"-^""'^""""({x};r,*)/a.(x„.,x,,i), (95) 

1=1 m=l 

n n 

Q^f--{{x}n,t) = J2Y1 5r'"^(x^,x,,o<ir"'"""'+'"""''({x};r)- (96) 



=1 m=l 



Here, if^i^a-""" ^ B"^ , and are defined as 



i?r""-""^({x}n,i) = 

.k=l 



(97) 



E(^Ax.+5Ax'Ji^n"^"^-""({X}„,t) 
k=l 

B'l'='{X,Y,t) = {Uf,{Y,t)wa{X,t)) - {Up{Y,t))) {wa{X,t)) , (98) 

2 

/«(Xi,X2,i) = ^(5;3,x, +9/3,x'JSf (Xi,X2,i). (99) 



k=l 
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It is worth noting that the assumption u G IKq imphes that i7„({X}„,t) = and therefore 
/„ (i)({X}„, t) = /„ (3)({X}„, t) = 0. However, we wih retain generahty and keep ah three terms. 
It should also be stressed that we are not assuming statistical independence between the velocity 
field na(x, and the velocity differences i(;a(x, x', t). On the contrary, we assume that the two are 
related to each other in the sense that the velocity field tiQ-(x, t) is forcing the velocity differences 
Wai^jx' ,t) via the sweeping interactions. However, the case of total statistical independence gives 
exactly /n({X}„,t) = /„ (i-)({X}„, i), so it is covered by our argument bellow. 

Let A be the scaling exponent of B"^^ such that 

Bf{RX,RY,t)r^(R/£o)\ (100) 

for scales in the inertial range. It immediately follows that la scales as 

Ia{RXi,RX2,t) ~ {R/£o)\l/R)giR). (101) 

Here, g{R) is a smooth function which represents departure from local homogeneity in the sense 
u € iKi. Without loss of generality, we associate the scaling exponent b to the function g{R). The 
contribution (1/-R) arises from the derivatives. Using a similar line of argument we see that Hn 
scales as 

//r""-""(«{X}n,i) ~ {R/iof-{l/R)UR), (102) 

where fniR) is also a smooth function representing departure from incremental homogeneity u G 
"Kq. We associate the scaling exponent to the function fn{R)- The three contributions to 
In{{X.}n,t) then scale as 

4,(l)(i?{X}„,t) ~ (i?/4)^"(l/i?)/n(i?), (103) 

4,(2)(iZ{X}„,t) ~ {R/io)'^-^{R/io)\l/R)giR), (104) 

4,(3)(i?{X}„,t) ~ {R/io)\R/eo)'^"''{l/R)fn-i{R), (105) 
and from power counting we find that the corresponding scaling exponents are 

An,i = Cn - 1 + ^n, A„,2 = Cn-1 + A - 1 + 6, and A„,3 = Cn-1 + A - 1 + a„_i. (106) 

Using the multifractal formulation, the contribution that supports the Holder exponent h gives 
(^n = nh + Z{h) , which gives the following evaluation for the scaling exponents 

An,l(/l) = (Cn - 1 + an) - (Cn+1 - I) = -h + a„, (107) 
A„,2(/i) = (Cn-1 + A - 1 + 6) - {Cn+1 - 1) = -2/i + A + /3, (108) 
An,3(/i) = (Cn-1 + A - 1 + a„_i) - (Cn+1 - I) = -2h + \ + a„_i. (109) 

Because the functions fniR) and g{R) are smooth, we can Taylor-expand them around R = 
and get, to first order, = 6 = 1. It is also reasonable to assume that A > since B'^^ involves a 
velocity difference. Prom these evaluations we find that the window for positive scaling exponents 
A„ is at least h G (0, 1/2). Admittedly, this is a rather narrow interval, even though it is sufficient for 
the downscale energy cascade of three-dimensional turbulence. However, the situation is probably 
a lot better than that. If we allow negative evaluations of R , which can be defined by reflecting 
the points {X.}n around their center of mass, we may expect that i? = is an extremum and 
therefore fniR) = g'{R) = 0. It is easy to show that, using the evaluation a„ = 6 = 2, we find 
that the window for positive scaling exponents A„ covers the entire range h G (0, 1) of local scaling 
exponents. Although this is somewhat encouraging, the real challenge is to determine what happens 
in reality and make a comparison of that against the speculative predictions given above. 
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6 Discussion and Conclusion 



In the original formulation of his theory, Kolmogorov assumed local homogeneity, local isotropy, 
and local stationarity in a non-Eulerian representation very similar to the quasi-Lagrangian rep- 
resentation of Bclinicher and L'vov. Frisch [7, 72] revised this argument by stating the same as- 
sumptions in the Eulerian representation. This is a decision that we agree with, because the energy 
spectrum and the structure functions are both Eulerian rather than Lagrangian quantities. Fur- 
thermore, as we have argued in the previous section, it is desirable to justify the elimination of 
the sweeping interactions theoretically, rather than hide the problem under a stronger definition 
of local homogeneity. Frisch [7] has also chosen to strengthen the assumption of self-similarity 
to make it possible to deduce all the scaling exponents (^n and obtain the prediction ^„ = n/3. 
Ultimately, this assumption needs to be replaced with a weaker assumption of self-similarity to 
permit intcrmittcncy corrections, and this is the approach followed in the L'vov-Procaccia theory 
[13,46,47,51-55]. Frisch himself proposed the multi- fractal hypothesis [7], which converges with 
the approach of Belinicher et al. in the papers [52-54] in a very interesting way. Finally, Frisch [7] 
made the very important observation that, in order to carry Kolmogorov's argument through, it is 
necessary to assume the existence of an anomalous energy sink. 

In the present paper we have shown that the assumptions of the Frisch framework are still not 
strong enough to prove the 4/5-law in the Eulerian representation. We have also shown that the 
problem of eliminating the sweeping interactions with a predictive argument remains open. The 
hypothesis u G ^i{A) can rectify both problems. Even better, starting from the more reasonable 
hypothesis u G %q{A), a rigorous proof that establishes A„ > would be sufficient to establish 
u € !Ki{A). A positive response to the question raised recently by Frisch [12] concerning the self- 
consistency of incremental homogeneity in the u G 5f o i-^) sense would also follow from the validity 
of the conjecture A„ > 0. This would be a fundamental breakthrough finally putting to rest the 
problem of the sweeping interactions that has concerned the community for the last 60 years. It 
would essentially establish that the sweeping interactions can be modeled as stochastic forcing 
acting only at large scales. Then we can simply drop from the Navier-Stokes equations the portion 
of the nonlinearity associated with the sweeping interactions, and build the entire statistical theory 
on the modified Navier-Stokes equations. We have also explained why it may not be desirable to 
use the Belinicher-L'vov quasi-Lagrangian formulation to go around the problem. The reason is 
that using the quasi-Lagrangian formulation requires the even stronger assumption u G !H* . 

We would also like to emphasize that our conclusions arc not a criticism of the L'vov-Procaccia 
theory based on the quasi-Lagrangian formulation [13,46,47,51-55]. In fact, as long as one's in- 
tention is to solve the problem of globally homogeneous turbulence, there is no issue whatsoever 
with respect to the sweeping interactions and the quasi-Lagrangian transformation, provided that 
the assumption of global homogeneity is many-time rather than one-time. On the other hand, it 
is desirable to move away from the assumptions of global homogeneity and global isotropy, which 
cannot be physically realized, and take steps towards building a theory based on the assump- 
tions of asymptotic incremental homogeneity u G J{q{A) and incremental isotropy in an Eulerian 
framework, as envisioned by Frisch. Our paper implies that the results of the L'vov-Procaccia the- 
ory [13,46,47,51-55] can be readily carried over and applied towards this goal, provided that the 
hypothesis A„ > is proved, and our proposal of section 5.2 rather than the quasi-Lagrangian 
transformation is employed. 
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A Quasi-Lagrangian representation of the Navier-Stokes equations 

We show how the quasi-Lagrangian transformation makes it possible to write the Navier-Stokes 
equations in terms of velocity differences, thereby eliminating the sweeping interactions. The reader 
should compare this argument with the derivation of the 4/5-law by Monin and Yaglom [73], to see 
that they are doing the same thing. 

We begin by noting that the Eulerian velocity field Wq(x, t) can be reconstructed from Vai'^o, io|x, t) 

as: 

ltQ,(x,t) = Uq,(x - p{xQ,to\t) + p{xo,to\t),t) = Va{xQ,tQ\x - p{xo , tQ\t) , t) . (A.l) 

To eliminate p{x.Q,to\t) we use 

Pa{xo,to\t)= dr Uaix + p{xo,tQ\T),T) = dr Va{xQ, tQ\x, t) , (A.2) 

Jto Jto 

and therefore, ^^(x, t) reads 

■Ua(x,t) =Va (xo,to\x- dx v(xo, to|x, t), . (A.3) 
Let RaixQ,to\x,t) be defined as 

i?a(xo,to|x,t) = (x)q - /" dr Va{y:o,to\x,T) = {x)a- Pa{y^O,to\t), (A.4) 

such that Ua{x,t) = Va{xQ,to\'R{xo,tQ\x,t),t). Also, define v^{xQ,to\t) = fQ(xo, to|xo, t) such that 
we may write concisely Wa = Va — v'^. It is easy to see that da,xRf3 = <^a/3 and dRa/dt = —v'^, 
and we may use these relations to show that the quasi-Lagrangian transformation preserves incom- 
pressibility, as follows: 

da,:>iWa = da^^{Va " V^^) = da,y.Va = {dp^^Va)^^^ = {di3,:,^Va){da,^Rf3) = da^^Ua = 0. (A.5) 

The key result is that the sweeping interactions are eliminated in the transformation of the 
material derivative itself. The show this, consider an arbitrary field U (x, t) and its quasi-Lagrangian 
transformation 'U,(xo, to|x, t) (where the fluid particle follows the Eulerian velocity field ). 
Prom the relation 

U{x, t) = U(xo, to|R(xo, io|x, t), t), (A.6) 
we find that 
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and 

da,^U = {dp,^U){da,^Rf3) = {dp,^U)Sap = {da,^U), (A.8) 

and it follows that: 

DU dU ^ dU ^ ^, 

-fZ- = ^ + Uada,xU = ^ + Wado,,^U. (A.9) 

Dt at ot 
This equation is identical to the unnumbered equation preceding equation (22.14) in Monin and 
Yaglom [73]. It is easy to see that since the material derivative is written in terms of velocity 
differences, if it is applied on ^^^(xo, to|x, we shall obtain an equation written exclusively in 
terms of velocity differences. 



B Evaluation of J[u\ 

In this appendix, we provide a detailed evaluation of the functional determinant J[u\ that we 
encounter in the derivation of the MSR theory for Lagrangian trajectories. The procedure was 
outlined briefly in L'vov and Procaccia [13]. However it is not as trivial as it seems. We have 
followed the outline and rederived the following more complete version of the proof: 

First, wc discretize time in Ai intervals and introduce the following notation: 



Pa = Pa{^0,to\tn) 
<'"^=««(xo+pS,tm)- 



(B.l) 



Each of the objects p^,u'^,Ua is a field that is a function of Xq only. Note that = 0. The 
governing equation for the Lagrangian trajectories field is equivalent to a set of the following 
discretized equations: 

n+l _ n 

^" ^" = <■ (B.2) 

There are, of course, many alternative discretizations to choose from. The rule is that, once we 
have chosen a discretization, we have to stay with it. We cannot switch to another scheme in the 
middle of the computations, for the sake of convenience. To evaluate J[u\ we proceed from the path 
integral definition: 



J[u]= [ l)pS[pa-Ql°g[p]up] 
= lim n / ^6 



neZ-{0} ' 
-oo 



n+l _ n 
ra ro 



At (B.3) 



limn^ rl 

n=l n=— oo 



where A„(At) and Bn{At) arc defined as 
B„{At) = jdpl5 



^n+l _ n 

At " 



(B.4) 
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Here, a is a normalization constant such that the product in (B.3) converges. Obviously, if such a 

constant exists, it will be unique. 

The An integral is easy to evaluate: 

An{At) = J dp^+^ AtS{pl+^ -pl- At<) = Ai. (B.5) 

To evaluate Bn we need to rewrite the discretized governing equation so that it is explicit with 
respect to p": 



(B.6) 



pl+^-pl = MUo.{^0 + Pl,tn) 

= At«„(xo + pl+^ - p^+^At + 0{Af),tn) 
= <+i'"Ai - (p^+i - p-p)dpul+''^At + 0{At'), 

therefore 

{5a(3 + dpu-^-^''^At){p^+' - p^p) = <+i'"At. (B.7) 
We proceed by employing the following change of variables: 

K = (^a/3 + 5/3<+''"A0p^. (B.8) 

The integral differentials arc transformed according to a determinant as follows: 

dRl = det((5a^ + 5^<+^'"A0dp3. (B.9) 

We will now show that incompressibility implies that the determinant is equal to 1. For brevity, 
introduce 

Map = Sap + dpul+^^^At. (B.IO) 

In the determinant expansion, every term other than M11M22M33 is O(At^) because it includes at 
least two off-diagonal factors each of which contributes a factor of At . It follows that: 

det M = M11M22M33 + 0{At^) 

= (1 + Atai + At92n^+''")(1 + Ata3tX3+''") + 0{Af) 

= 1 + At(ai<+^'" + 92%^''" + 93%+'''') + 0{At^) 
= l + 0{At^). 

Note that in the last step wc employed the incompressibility condition. It follows that dp"' = 
dR" . We may now proceed and evaluate the integral B^- 



Bn 



{At) = I dp"AtS[pl+'-pl-Atul] 

= Atj dp^^ 5{{5o.p + a^<+i'"At)(p^+i - pj) - ul+^'^At) 

= Atil + 0{Ae)) ^^-^^^ 

X j dR^ 5{{5ap + dpul+^^^At)pl+' -Rl- Ai<+i'-) 
= Ai + 0(Af^). 
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In the last step, the crucial requirement is that and n"^ should not depend on and 

therefore If we set the normalization constant a = At, then J[u\ evaluates as: 



n=l 



Note that O(At^) contributions to the integrals A„(Ai) and Bn{A.t), which we have disregarded, 
would vanish anyway after taking the limit At — ^ 0, so they can be safely ignored with impunity. 



C Sweeping interactions under a gaussian mean field 



We exploit the following mathematical result: if /a(xi,ti) is a Gaussian stochastic field, the 
ensemble averages of the form (/q(xi, ti)i?[/]) can be evaluated for any analytic functional R[f] by 
the following integral 

(/a(xi,ti)it;[/]) = (/«(xi,ii))(i?[/])+ j d^2dt2 {fa{^l,h)ffs{^2,t2))c{j^^) ^ (d) 



where 



(/a(xi,tl)//3(x2,t2))c = (/a(xi,tl)//3(x2,t2)) - (/a(xi, tl)) (//3(x2, ^2)) • 



(C.2) 



This a generalization of Gaussian integration by parts, a technique attributed by Prisch [7] to 
Novikov [31], Donsker [106] and Furutsu [107]. 

We begin the proof by defining the following correlation functions: 



Uap{^l,ti;y.2,t2) = («a(xi,il)u/3(x2,t2)) - («a(xi,tl)) {up{yL2M)) 



5a...a„/J(|x}„,Y,t) = ((lX;3(Y,t)) 



.1=1 



,fe=l 



and also the following response functions: 



Rr 



_ I dWa{X,ti) 

\Suf3{y,t2) 



S 



.1=1 



(C.3) 

(C.4) 
(C.5) 

(C.6) 
(C.7) 



Here, we disregard the fact that tia(x, t) and Wq,(x, x',t) are related by definition and assume that 
the only effect of Ua{x,t) on Wa{x,x' ,t) is via the sweeping interaction. We also assume that the 
velocity field Ua{x,t) is delta-correlated which implies that 



^^a/3(xi,ti;X2,i2) = C/a/3(xi,X2)5(il -t2). 

We begin by splitting /„({X}„,t) into two terms 

-^n({X}n,t) = /„^(i)({X}„,t) + /n,(2+3)({X}n,i)) 



(C.8) 

(C.9) 
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with 7„_(i)({X}„,t) given by 



k=l 

= {U(,i{xu,t)) 

Here we have used the incompressibihty condition. 



]Ju;a,(X/,t) 



./=i 



.fc=i 



(C.IO) 

(C.ll) 
(C.12) 



fc=i 

The remaining contribution to /„({X}„,t) reads 

^ n n 
fe=l i=l 

Using Gaussian integration by parts we may write 

Bf{X, Y,t) = J dzdT RayiX, t; z, r)[C/^^(y, t; z, r) + C/^^tCy', ^)] 
= y dz i?a7(X,i;z,i)[?7/3^(y,z) + C/^^(y',z)], 

and 



(C.13) 



(C.14) 



(C.15) 
(C.16) 



^•-"^({X}„,Y,t)= j dzdri?-i-"^({X}„,t;z,r)[C/^^(y,t;z,r) + C/^^(y',t;z,r)] (C.17) 
= y"dzi?^i-"»^({X},,i;z,t)[[/^^(y,z) + C/^^(y',z)]. (C.18) 

The key step is to note that 



K--""'^({X}„,i,y,i) 



n 



k=l 



E F^i-"''-^"''--"" ({x}^)i?.,Mx., y, t). 



(C.19) 
(C.20) 
(C.21) 



k=l 



Here we exploit the fact, first pointed out in Ref. [47], that the variational derivative {5wa^. {^k,'t))/i^Ui3{y, t)) 
is not correlated with the velocity differences w;^, (X; , t) because no time is being allowed for the 
interaction to develop a correlation. This relationship between the response functions implies a 
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corresponding relationship between B„({X}„, Y,t) and Bi{X.k,Y ,t): 



.k=l 



[Up^{y,z) + Up^{y',z)] 
(C.22) 

(C.23) 



k=l 



= E^n-r'"'"'+'''""({X}^)^r^(Xfc,Y,t). (C.24) 
k=l 

It immediately follows that 

^ n n n 

Cm" ^^^>- = 2^ E E E (^/3,x. + %x', )F:rr'"-^"'"+"""" ({x}-)sr^(x„, x,, 



A;=l 1=1 m=l 



which can be broken down to 

n n 



C(2f""({x}n,t) = E E i^'^r""-^""'+"-""({x};r)^.^(x„»,x,,i), 

;=1 m=l 

n n 

7;;j--"({x}„,t) = E E sr'(x™,x,,tXir-^"-+"-""^({x}-), 



(C.25) 

(C.26) 
(C.27) 



=1 m=l 



with 



/a(Xi,X2,i) = E(9Ax, +5/3,x'JSf (Xi,X2,i). 



(C.28) 



k=l 
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